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Abstract. The results of single-dish observations of low- and high- J transitions of selected molecules from pro- 

>^ , toplanetary disks around two T Tauri stars (LkCa 15 and TW Hya) and two Herbig Ae stars (HD 163296 and 

^ ■ MWC 480) are reported. Simple molecules such as CO, ^^CO, HCO+, CN and HCN are detected. Several lines 

of H2CO are found toward the T Tauri star LkCa 15 but not in other objects. No CH3OH has been detected 

down to abundances of 10~^ — 10~* with respect to H2. SO and CS lines have been searched for without suc- 

^— ^ ' cess. Line ratios indicate that the molecular emission arises from dense (10^-10* cm~^) and moderately warm 

^^ I (r ^ 20-40 K) intermediate height regions of the disk atmosphere between the midplane and the upper layer, 

f^ . in accordance with predictions from models of the chemistry in disks. The sizes of the disks were estimated from 

model fits to the ^^CO 3-2 line profiles. The abundances of most species are lower than in the envelope around 

the solar-mass protostar IRAS 16293-2422. Freeze-out in the cold midplane and photodissociation by stellar and 

interstellar ultraviolet photons in the upper layers are likely causes of the depletion. CN is strongly detected in all 

O ' disks, and the CN/HCN abundance ratio toward the Herbig Ae stars is even higher than that found in galactic 

photon-dominated regions, testifying to the importance of photodissociation by radiation from the central object 

C/j ■ in the upper layers. DCO^ is detected toward TW Hya, but not in other objects. The high inferred DCO^/HCO^ 

' ratio of ~0.035 is consistent with models of the deuterium fractionation in disks which include strong depletion 

of CO. The inferred ionization fraction in the intermediate height regions as deduced from HCO"*" is at least 

10^^^ — 10"^'', comparable to that derived for the midplane from recent H2D''' observations. Comparison with 

the abundances found in cometary comae is made. 
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1. Introduction of a disk, and after the collapse and accretion onto the 

star cease, planets and comets can form by accumulat- 

The protoplanetary disk phase constitutes a key period -^^ ^^^^^^^ ^^^ ^^^-^ material on timescales of a few 

in the evolution of matter between the young protostellar ^.^^.^^ ^^^^^^ ^^ ^^^ Lissauer HHSI Beckwith & Sargent 

and the mature planetary system stages. Before their in- ^^^ Wuchterl et aLEnnOJ. Surveys from the near-infrared 

corporation mto comets and large bodies, the gas and dust ^^ ^^^ millimeter wavelength range have shown that a 

could have participated in a complex chemistry within the 1 r _,■ rim -n- u o i-i 

. ^ ^ ... . large traction 01 1-10 million year old bun-like pre-mam- 

disk. Studies of the chemistry in disks are therefore impor- , 11 j-i-t.^i- ^i- / 

■' ^ sequence stars harbors a disk m Kepierian rotation (e.g., 

tant to quantify the chemical composition of protoplane- g^^kwith et al. [IS9- The masses of these disks (0.001 - 

■^ ■ 0.1 Mq) is sufficient to form a few giant gaseous planets. 

The chemical composition of the envelopes around young Single-dish and interferometric observations of molec- 

protostars is now known with increasing detail thanks to ^^^r species other than CO are starting to reveal the 

the combination of rapid advances in detectors and an- chemistry m disks around classical T Tauri s_tar^(Dutrey, 

tennatechnologyandimprovedmodels(e.g.,vanDishocck GuiUoteau & Guclin 1997: Kastner et al. 1997; Simon, 

& Blake in^, Langer et al.^m- Part of this gas and D^^t^^^y ^ GuiUoteau ,2000. Duvert etaL ,20001 van 

dust settles around the pre-main-sequence star in the form Zadelhoff et al. 200]UAikawa et al. EOmiQi et al. 

EOOa Dartois et al. HUM Kessler et al. HUM Wilner 

Send offprint requests to: Ewine F. van Dishoeck, e-mail: et al. .2003,) . The low- J rotational transitions of simple 

ewineOstrw. leidenuniv.nl molecules (HCN, CN, HNC, H2CO, HCO+, CS, ...) are 
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detected, but their abundances relative to H2 are inferred 
to be orders of magnitude lower than those observed in 
dark clouds. The prevailing explanation of this depletion 
involves a combination of freeze-out of the molecules on 
grain surfaces in the cold midplane and their photodis- 
sociation by ultraviolet and/or X-rays in the upper atmo- 
sphere of disks (see Aikawa et al.[T?)99a 2002; Bergin et al. 
I2()()H|I . The abundances are enhanced in the intermediate 
height regions, which are warm enough for the molecules 
to remain in the gas phase. Photodesorption induced by 
ultraviolet radiation (Willacy & Langer I^OOOI Westley et 
nl. HM^ or X-rays (Najita, Bergin & UUom lMTHl can fur- 
ther populate the upper layers with molecules evaporated 
from dust grains. 

We present here the results of a survey of several low- 
and high-J molecular transitions observed toward two 
classical T Tauri stars (LkCa 15 and TW Hya) and two 
Herbig Ae stars (MWC 480 and HD 163296) using single- 
dish telescopes. In particular, organic molecules such as 
H2CO, CH3OH and HCN and deuterated species were 
searched for. The comparison of the two types of objects 
allows the influence of the color temperature of the radia- 
tion field on the chemistry to be studied. There are several 
advantages in observing high-J transitions over the lower- 
J ones. First, detections of CO J = 6^5 and H2 show the 
presence of a warm upper surface layer in protoplanetary 
disks whose temperature is higher than the freeze-out tem- 
perature of most volatile molecules (van Zadelhoff et al. 
IMm Thi et al. l2()()l|l . Combined with the high densities in 
disks, this allows the mid- J levels to be readily populated. 
Models of flaring disks predict that the upper layer fac- 
ing directly the radiation from the central star can extend 
out to large radii (Chiang & Goldreich ll997l D'Alessio et 
al. ll999|l . Second, by observing at higher frequencies with 
single dish telescopes, the lines suffer less beam dilution 
entailed by the small angular size of disks, typically 1-3" 
in radius, than at lower frequencies. Also, confusion with 
any surrounding low-density cloud material is minimized. 

The results for the different molecules are compared 
to those found for protostellar objects, in particular the 
solar-mass protostar IRAS 16293-2422, which is consid- 
ered representative of the initial cloud from which the Sun 
and the solar nebula were formed. This so-called Class 
object (Andre et al. I2[)()[)|l is younger than the proto- 
planetary disks studied here, only a few x 10^ yr, and its 
chemistry is particularly rich as shown by the number of 
species found in surveys in the (sub)millimeter range (e.g., 
van Dishoeck et al. 119951 Ceccarelli et al. 120011 Schoier et 
al. 120021 Cazaux et al. 120031 and references therein). The 
similarities and differences in the chemical composition 
between IRAS 16293-2422 and the protoplanetary disks 
can be used to constrain the chemical models of disks. 

At the other extreme, the results for disks can be com- 
pared with those found for objects in our solar system, in 
particular comets. This will provide more insight into the 
evolution of matter from the protoplanetary disk phase to 
planetary systems. Unfortunately, the chemical composi- 
tion of the large bodies in our solar system has changed 



since their formation 4.6 Gyr ago. For example, solar radi- 
ation triggers photochemical reactions in the atmospheres 
of planets, and the release of energy from the radioactive 
decay of short-lived elements such as ^^Al causes solids 
to melt. Comets, however, could have kept a record of the 
chemical composition of the primitive solar nebula because 
they spent much of their time in the cold outer region of 
the Solar System (the Oort cloud) since their formation 
(Irvine et al. 120001 Stern '20JJ3. Comparison of cometary 
D/H ratio and the CH3OII abundances with those in disks 
are particularly interesting. 

This paper is organized as follows. In § 2, the character- 
istics of the observed objects are summarized. In § 3, the 
observational details are provided. The results are given 
in § 4 where a simple local thermodynamical equilibrium 
(LTE) and statistical equilibrium analysis is performed. In 
this section, we also derive several disk characteristics by 
fitting the ^^CO 3-2 lines. In § 5, the molecular abundance 
ratios are discussed. In particular, the CN/HCN ratio can 
trace the photochemistry whereas the CO/HCO'^ ratio is 
a tracer of the fractional ionization. Finally, a discussion 
on the D/H ratio in the disks compared with that found 
in comets or other star-forming regions is presented (see 
also van Dishoeck et al. l20()3|l . 



2. Objects 

The sources were selected to have strong CO J = 3 ^ 2 
fluxes and the highest number of molecular lines detected 
in previous observations (Qi 120011 Thi et al. 120011 van 
Zadelhoff et al. "20Q1I). LkCa 15 is a solar mass T Tauri 



star located in the outer regions of the Taurus cloud. Its 
age is estimated to be ~10 million years, although Simon 
et al. (2000) argue for an age of only 3-5 million years. 
LkCa 15 is surrounded by a disk whose mass is estimated 
to be around 0.03 Mq, although a higher mass has been 
obtained from the fitting of its spectral energy distribution 
(SED) (Chiang et al. 2001) . LkCa 15 is one of the strongest 
millimeter emitting sources in the sample of T Tauri stars 
surveyed by Beckwith et al. H1990|l along with GG Tau 
and DM Tau. 

TW Hya forms part of a young association of stars 
that has been discovered only recently and is located at 
only ~56 pc (Webb et al. lTO^ . TW Hya itself is a clas- 
sical isolated T Tauri star with a high X-ray flux and a 
large lithium abundance. Its large Ha equivalent width 
is indicative of active disk accretion at a rate of ~ 10~^ 
Mq yr^^ (Kastner et al. l2002|l . Despite its relatively high 
age (~15 Myr), TW Hya is surrounded by a disk of mass 



3x10-2 M0 (Wilner et al. 



seen nearly face-on 



(Weintraub et al.HnHi Krist et al lMIUI Zuckerman et al. 

mm- 

MWC 480 and HD 163296 were chosen to be repre- 
sentative of Herbig Ae stars. These two objects have the 
strongest millimeter continuum emission, with disk masses 
similar to those around the two T Tauri stars. All se- 
lected objects show gas in Keplerian rotation as revealed 
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Table 1. Stellar characteristics" 



Star 



SpT 



(J2000) 



d 
(J2000) 



D 

(pc) 



M. 
(Mq) 



Log(i./LQ) 



LkCa 15 
TW Hya 
HD 163296 

MWC 480 



K7 
K8Ve 
A3Ve 
A3ep+sh 



04 39 17.8 
11 01 51.91 
17 56 21.26 
04 58 46.27 



+22 21 03 
-34 42 17.0 
-21 57 19.5 
+29 50 37.0 



140 
56 
122 
131 



0.8 
1.0 
2.4 
2.2 



-0.27 
-0.60 
+1.41 
+1.51 



Age 
(Myr) 



11.7 
7-15 
6.0 
4.6 



° See Thi et al. (2001) for references. The ages are highly uncertain (see also Simon et al. l2000ll . 



Table 2. Disk characteristics 



Star 


Disk Mass" 


Radius 


Diameter 


Inclination 




Rcf. 




(W-'Mq) 


(AU) 


(") 


n 






LkCa 15 


3.3 ± 1.5 


425 


6.2 


57+10 


1 




TWHya 


3.0 ± 2.0 


200 


7.0 


<10 


2, 


3, 4, 5, 6 


HD 163296 


6.5 ± 2.9 


310 


5.0 


32+5 


7 




MWC 480 


2.2 ± 1.0 


695 


10.4 


30+5 


7 





" Total gas + dust mass computed from millimeter continuum flux using a dust opacity coefficient KA=0.01(1.3mm/A) cm^ 
g~^ and assuming a standard gas-to-dust mass ratio of 100 

Ref erence s. (1) Qi et al. (PTM)l; (2 ) Weintraub et al. C^l; (3) Krist et al. lUEl; (4) Weinberger et al. I^FM^ : (5) Wilner et 
al. lUTOni : (6) Calvet et al. J^HH^ : (7) Mannings et al. (tB^ 



by CO interferometric observations (Qi l200ll Mannings & 
Sargent 1997) ■ 

The stellar characteristics of the four objects are given 
in Tableland the disk characteristics in Tabled In this 
paper, inclination is defined such that 0° is a disk viewed 
pole-on and 90° edge-on. Detailed modeling of the SED 
of these objects suggests that the disks are in a state of 
dust settling, especially LkCa 15 (Chiang et al. 12001(1 . The 
objects were chosen to be isolated from any cloud material 
to avoid contamination of the single-dish data. 



3. Observations 

The observations were performed between 1998 and 2000 
with the James Clerk Maxwell Telescope (JCMT)i lo- 
cated on Mauna Kea for the high-J transitions (850 /im 
window) and with the 30-m telescope of the Institut de 
Radioastronomie Millimetrique (IRAM) at Pico Veleta for 
the lower J lines (1 to 3 mm). At both telescopes, the ob- 
servations were acquired in the beam-switching mode with 
a throw of 120" at IRAM and 180" at the JCMT in the 
azimuth direction. The observations suffer from beam di- 
lution owing to the small projected sizes of the disks of 
the order of 5-10" in diameter, compared to the beam 
size of the JCMT (14" at 330 GHz) and IRAM (11.3" at 
220 GHz). The data were reduced and analyzed with the 
SPECX, CLASS and in-house data reduction packages. 



^ The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre in Hilo, Hawaii on behalf of the 
Particle Physics and Astronomy Research Council in the 
United Kingdom, the National Research Council of Canada 
and The Netherlands Organization for Scientific Research. 



The JCMT observations made use of the dual polariza- 
tion B3 receiver (315-373 GHz) and were obtained mostly 
in November-December 1999. The antenna temperatures 
were converted to main-beam temperatures using a beam 
efficiency of ?7,nb=0.62, which was calibrated from obser- 
vations of planets obtained by the staff at the telescope. 

The data were obtained in single sideband mode with 
the image side band lines reduced in intensity by about 
13 dB (i.e. by a factor of ~ 20). The sidebands were cho- 
sen to minimize the system temperature and to avoid any 
unwanted emission in the other sideband. The integration 
times range from 5 minutes for the bright ^^CO J = 3— >2 
lines to 8 hours for the faint lines to reach a r.m.s. noise 
STthcr of 10-20 mK after binning. The backend was the 
Digital Autocorrelator Spectrometer (DAS) set at a reso- 
lution of-- 0.15-0.27 km s'^ (see TablesElandEj), and sub- 
sequently Hanning-smoothed to 0.3-0.6 km s^^ in spectra 
where the signal-to-noise ratio is low. Pointing accuracy 
and focus were regularly checked by observing planets, 
and was found to be accurate to better than 3" rms at the 
JCMT. 

The estimated total r.m.s. error 6T at the JCMT 
associated with each line is given by the relation (e.g., 
Papadopoulos & Seaauist 11^98) : 



ST 

Tmh 



ST 

Tmh 



thcr 



ST 

Tmh 



drift 



ST 
Tnih 



syst 



1/2 



(1) 

where the first term on the right-hand side of the relation 
expresses the ratio between the thermal r.m.s. tempera- 
ture and the main-beam peak temperature averaged over 
A^ch channels and with a baseline derived from A'^bas chan- 
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nels: 



5T 

Tmh 



ST, 



thcr 



the 



A 



Nu 



N, 



ch 



A^bas^ch 



1/2 



(2) 



where iJTthcr is the thermal noise per channel and T^ is the 
antenna temperature given in Tables |31 and 01 The beam- 
switching method gives extremely flat baselines such that 
A^bas >> Nch and 



ST 

Tmh 



ST 



tlicr 



thcr 



T* 



(3) 



The data are subsequently binned to A'ch = 2 or 4, and a 
line detection is claimed whenever 



ST 

Tmh 



<0.3 



(4) 



the 



The second therm in Eq. Q] concerns the drift of 
the temperature scale due to all errors of stochastic na- 
ture, and therefore includes any temperature variation of 
the cold loads or any fluctuation of atmospheric opac- 
ity. Measurements of spectral standard sources just be- 
fore or after the source observations allow an estimate of 
this drift, which is generally found to be 10-15% and may 
be up to 20-25% in difficult parts of the atmospheric win- 
dow (e.g. 1120+ amd N2H+ lines) depending on the condi- 
tions. The lines for which calibration sources are available 
from measurements by the JCMT staff are mentioned in 
Table Ef. 

The last term encompasses the systematic error, whose 
main contributors are the uncertainty in the value of the 
beam-efficiency and pointing errors. As noted above, the 
pointing at the JCMT was found to be accurate to bet- 
ter than 3" at 345 GHz. Differences in beam efhciencies 
and pointing should also be reflected in the spectral stan- 
dard observations, which generally agree within 10-15% as 
noted above. This last term is therefore estimated to con- 
tribute at most 10-20%. Taking into account all possible 
sources of errors, the overall calibration uncertainties can 
be as high as 30-40% for a line detected with 3a in a dif- 
ficult part of the spectrum, whereas it is of order 20-25% 
for high S/N lines for which spectral standards have been 
observed. 

As discussed by van Zadelhoff et al. H2001|l . our HCO+ 
J = 4—3 intensity measured in 1999 is a factor of 3 weaker 
than that obtained by Kastner et al. i|1997|l . More recently, 
we have re-observed the HCO"*" line in May 2004 and find 
intensities on two days which agree with those of Kastner 
et al. within 10-20%. For comparison, ^^CO 3-2 spectra 
taken in 1999, 2000 and 2004 are consistent within 10% 
with the Kastner et al. results taken in 1995 with a differ- 
ent receiver, as are the HCN 4-3 and CN 3-2 results. Thus, 
only the 1999 HCO"^ result appears anomalously low, per- 
haps due to unusually large pointing errors during those 



observations related to the JCMT "tracking error" prob- 
lem '^, unless the ion abundance is variable. We use only 
the new 2004 data in our analysis. Note that the H^'^CO+ 
and DCO"'' data were taken only 1 week apart so that the 
analysis of the DCO+/HCO+ ratio should not suffer from 
any potential long-term variability. Further monitoring of 
the HCO+ line is warranted. 

The IRAM-30m observations were carried out in 
December 1998 using the 1-3 mm receivers. The weather 
conditions were excellent. The three receivers and a 
splitable correlator were used to observe simultaneously 
lines at 1.3, 2 and 3 mm. The receivers were tuned single- 
sideband. Image band rejection was of the order of 10 dB. 
Forward (Foff ) efficiencies were measured at the beginning 
of each run and have been found to be consistent with 
standard values. We measured Fcff=0.9, 0.82 and 0.84 at 
100, 150 and 230 GHz respectively. The derived beam ef- 
ficiencies (?7mb = Bcff/Fos) are 0.57, 0.69, and 0.69 at 1, 2 
and 3 mm respectively using main-beam efficiencies (-Boff) 
provided by the IRAM staff. The pointing and focusing 
accuracy were regularly checked to ensure pointing errors 
< 3" (r.m.s.) by observing planets and quasars. TW Hya 
is unfortunately located too far south to observe with the 
IRAM 30m telescope. 



4. Results 

4.1. General characteristics 

The measured antenna temperatures and thermal noise 
per channel width are summarized in Tables O and 01 The 
spectra are displayed in Fig.[21to01on the main-beam tem- 
perature scale for the four sources. ^^CO J = 3^ 2 and 
^^CO J = 3 ^ 2 are detected toward all objects. Apart 
from TW Hya, the profiles of the ^^CQ J = 3 ^ 2 spec- 
tra are double-peaked with peak separations of ^ 2 km 
s"^ for both three objects. The ^^CO J = 3 -> 2 spec- 
trum of MWC 480 shows a profile with slightly different 
peak strengths. However, the level of asymmetry is not 
signifcant compared to the noise. ^^CO J = 3^-2 obser- 
vations obtained with 30"offsets and position-switching 
to an emission-free position are shown in Fig. \S[ for the 
four objects. The maps around LkCa 15, TW Hya and 
MWC 480 confirm that these objects are isolated from 
cloud material. The observations at offset positions from 
HD 163296 show emission at velocities shifted compared 
with the velocity of the star. The extinction to HD 163296 
is sufficiently low that the extended low density emission is 
unlikely to arise from a foreground cloud. The offset emis- 
sion is only seen in ^■^CO, not in ^^CO or other molecules. 
Lines arising from high- J transitions require high critical 
densities and are therefore not likely to come from a low 
density cloud. A more complete discussion on the possible 
contamination by foreground and/or background clouds is 
given in Thi et al. (;2001|) . 



see also 'http : //www . j ach . hawai i . edu/ JACpublic/ JCMT/ 1 'http : //www . j ach . hawaii . edu/ JCMT/Facility_descript ioii/| 

Heterodyne.observing/Standards/ Pointing/problem.transit .html 
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Table 3. Antenna temperatures and noise per channel width for the T Tauri stars 





Telescope 


dv" 


Tl 


(^rthcr 






(km s-i) 


(K) 


(K) 




LkCa 


15 






'^CO J = 2->1 


IRAM30m 


0.10 


0.402 


0.159 


'^CO J = 3->2 


JCMT 


0.13 


0.258 


0.051 


"CO J = 3^2 


JCMT 


0.28 


0.074 


0.026 


C'*0 J = 2^1 


JCMT 


0.21 


<0.084 (2cr) 


0.042 


C'*0 J = 3^2 


JCMT 


0.14 


<0.054 (2cr) 


0.027 


HCO+ J = 4^3 


JCMT 


0.26 


0.055 


0.018 


H"CO+ J = 4->3 


JCMT 


0.13 


<0.078 (2cr) 


0.039 


DCO+ J = 5-*4 


JCMT 


0.13 


<0.048 (2cr) 


0.024 


CN J = 3|^2| 


JCMT 


0.27 


0.070 


0.018 


HON J = 4->3 


JCMT 


0.26 


0.051 


0.022 


CS J = 7-^6 


JCMT 


0.27 


<0.038 (2cr) 


0.019 


H2CO J = 2l2^1ll 


IRAM30m 


0.16 


0.021 


0.008 


Ha CO J = 3o3^2o2 


IRAM30m 


0.10 


0.025 


0.021 


H2CO J = 322 ^221 


IRAM30m 


0.10 


<0.036 (2cr) 


0.018 


H2CO J = 3i2^2ii 


IRAM30m 


0.10 


0.066 


0.016 


H2CO J = 5l5-^4l4 


JCMT 


0.27 


0.030 


0.016 


CH3OH J = 2x->1k 


IRAM30m 


0.12 


<0.014 {2a) 


0.007 


CH3OH J = 42^3i E+ 


IRAM30m 


0.11 


<0.036 {2a) 


0.018 


CH3OH J = 5k->4k 


IRAM30m 


0.09 


<0.036 (2a) 


0.018 


N2H+ J = 4^3 


JCMT 


0.12 


<0.454 (2cr) 


0.227 


H2D+ J = li„^lii 


JCMT 


0.12 


<0.372 (2cr) 


0.186 


TWHya 


^^CO J = 3^2'' 


JCMT 


0.13 


1.853 


0.066 


i^CO J = 3->2" 


JCMT 


0.13 


1.679 


0.241 


"CO J = 3^2 


JCMT 


0.14 


0.185 


0.046 


HCO+ J = 4^3'' 


JCMT 


0.13 


0.656 


0.230 


HCO+ J = 4->3" 


JCMT 


0.13 


1.197 


0.087 


H"CO+ J = 4^3 


JCMT 


0.13 


0.050 


0.016 


DCO+ J = 5^4 


JCMT 


0.13 


0.104 


0.023 


CN J = 3|^2| 


JCMT 


0.14 


0.597 


0.075 


HCN J = 4^3 


JCMT 


0.13 


0.369 


0.080 


H"CN J = 4^3 


JCMT 


0.13 


<0.056 (2cr) 


0.028 


HNC J = 4^3 


JCMT 


0.13 


<0.089 {2a) 


0.049 


DCN J = 5-^4 


JCMT 


0.13 


<0.068 {2a) 


0.034 


H2CO J = 3i2-^2ii 


JCMT 


0.10 


<0.088 {2a) 


0.044 


H2CO J = 5l5^4l4 


JCMT 


0.13 


<0.074 {2a) 


0.037 


CH3OH J^Ik-'Qk 


JCMT 


0.14 


<0.046 {2a) 


0.023 


N2H+ J = 4->3 


JCMT 


0.12 


<0.710 (2a) 


0.355 


H2D+ J^lio^lll 


JCMT 


0.12 


<0.596 (2a) 


0.298 


SO J = 88^77 


JCMT 


0.13 


<0.202 (2a) 


0.101 



° dv is the spectral resolution of the data in km s ^ ; TJ is the antenna temperature (not corrected for beam efhciency) and 
(Either is the thermal noise per dv channel. ^ Data taken in November 1999. '^ Data taken in May 2004. 



High- J lines of various molecules are detected in 
the disks. Lines with high signal-to-noise ratio toward 
LkCa 15, HD 163296 and MWC 480 show a double- 
peak structure corresponding to emission from a disk in 
Keplerian rotation viewed under an inclination angle i 
(Beckwith & Sargent 1993). The line profiles observed to- 
ward TW Hya are well fitted by a single gaussian, con- 
sistent with a disk seen almost face-on. The profiles show 
no evidence of extended velocity wings characteristic of 
molecular outflows in any of the objects. The velocity in- 



tegrated main-beam temperatures for the four sources are 
summarized in Table El This table includes the energy of 
the upper level of the transitions, their critical densities 
and frequencies, the telescope at which they were observed 
and the beam size. The critical densities are defined as 
Her = j4ui/X]/ Quij where Aui is the Einstein A coefficient 
of the transition u — > I and qui the downward rate co- 
efficient. They have been computed in the optically thin 
limit at 100 K using the molecular data listed in Jansen et 
al. (1994) and Jansen (1995). For optically thick lines, the 
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Table 4. Same than Table Ol but for the Herbig Ae stars. 





Telescope 


dv 


n 


-^thor 






(kms~') 


(K) 


(K) 


HD 163296 


'''CO J = 3^2 


JCMT 


0.23 


0.920 


O.ISO 


"CO J = 3^2 


JCMT 


0.28 


0.225 


0.052 


HCO+ J = 4->3 


JCMT 


0.26 


0.128 


O.OSO 


CN J = 3|^2| 


JCMT 


0.27 


0.115 


0.018 


HCN J = 4^3 


JCMT 


0.26 


<0.074 {2a) 


0.037 


H2CO J = 5l5^4i4 


JCMT 


0.27 


<0.204 (2a) 


0.102 


H2CO J = 2l2-»lll 


IRAMSOm 


0.17 


<0.018 {2a) 


0.009 


H2CO J = 3o3~>2o2 


IRAMSOm 


0.11 


<0.030 {2a) 


0.015 


H2CO J = 3i2->2ii 


IRAMSOm 


0.10 


<0.018 {2a) 


0.009 


CH3OH J = 2k^1k 


IRAMSOm 


0.12 


<0.014 {2a) 


0.007 


MWC 480 


i^CO J = 3^2 


JCMT 


0.13 


0.498 


0.073 


"CO J = 3-^2 


JCMT 


0.28 


0.102 


0.016 


C'^O J = 3->2 


JCMT 


0.14 


<0.062 (2a) 


0.031 


HCO+ J = 4^3 


JCMT 


0.26 


0.052 


0.017 


DCO+ J = 5^4 


JCMT 


0.13 


<0.106 {2a) 


0.053 


CN J = 3|-^2| 


JCMT 


0.27 


0.038 


0.017 


CS J = 7^6 


JCMT 


0.27 


<0.040 {2a) 


0.020 


HCN J = 4->3 


JCMT 


0.53 


<0.024 {2a) 


0.012 


H2CO J = 3i2^2ii 


JCMT 


0.43 


<0.030 {2a) 


0.015 


H2CO J = 5l5->4i4 


JCMT 


0.27 


<0.032 {2a) 


0.016 


H2D+ J=ll0^1ll 


JCMT 


0.12 


<0.090 {2a) 


0.045 


N2H+ J = 4^3 


JCMT 


0.12 


<0.126 {2a) 


0.063 


H2CO J = 2i2^1ll 


IRAMSOm 


0.21 


<0.016 (2a) 


0.008 


H2CO J = 322^221 


IRAMSOm 


1.37 


<0.006 (2a) 


0.003 


CH3OH J = 42-»3i E+ 


IRAMSOm 


1.37 


<0.006 (2a) 


0.003 



critical densities are lowered by roughly the optical depth 
of the line. 

The upper limits are computed assuming a main beam 
temperature which corresponds to twice the r.m.s. noise 
level in a 0.3 km s^^ bin and a line profile similar to that 
derived from fitting the ^'^CO J = 3— »2 lines. 

The ion HCO"'" is detected in all sources. Toward 
TW Hya, H^'^CO+ is also seen, and the ratio of integrated 
fluxes HC0+/H13C0+ of 24 is lower than the interstellar 
isotopic ratio [i2c]/[i3C] ~ 60, indicating that HCO+ hue 
is optically thick (van Zadelhoff et al. l2()()l|l . Emission from 
the CN radical in the J = 3 — > 2 line is surprisingly strong 
in all objects compared to HCN J = 4^3. HNC J = 4-^3 
was only searched toward TW Hya and has not been de- 
tected. The line intensities for some species in this object 
differ with previous observations by Kastner et al. (1997) 
(see van Zadelhoff et al. 2001, §3). In general, the two 
Herbig Ae stars display a less rich chemistry than the two 
classical T Tauri stars. In particular, HCN is not detected 
in either source in our observations. Qi H2001|l . however, 
reports detection of HCN J = 1 ^ toward MWC 480 
with the Owens Valley Millimeter Array (OVRO). 

Several lines of H2CO are seen toward LkCa 15 with 
the IRAM 30-m and JCMT, but not toward the other 
three disks. Deep searches for various CH3OH lines with 
the IRAM-30m and JCMT down to very low noise levels 



did not yield any detections. So far, CH3OH has only been 
seen toward LkCa 15 through its J = 5^^ 4j^ and J = 
423^ 3i3 lines using OVRO (Qi OUUTIl . This illustrates 
the power of interferometers to detect minor species in 
disks, because, with a beam diameter of l"-2", the fluxes 
are much less beam-diluted, in addition to a potentially 
larger collecting area when a great number of dishes is 
available. 

No CS J — 7^6 line nor lines of SO2, some of which 
occur fortuitously in other settings (e.g., near H2CO 351 
GHz), were detected toward LkCa 15. A deep limit on SO 
was obtained toward TW Hya (Fig. |21 and Table El . 

Finally, DCO"*" is detected for the first time in a disk, 
as reported by van Dishoeck et al. H2003|l . The DCO+ 
J = 5 — > 4 line is observed toward TW Hya with a strength 
similar to that of H^^CO+ J = 4^-3, but the line is not 
detected toward LkCa 15 and MWC 480, where Hi3C0+ 
is also not seen. There is a hint of a feature near the DCN 
J = 5 ^ 4 line toward T W Hya, but it is offset by a few km 
s~^. Deeper integrations or interferometer data are needed 
to confirm this. Note also the high critical density needed 
to excite the DCN J = 5^4 line (~ 4.8 lO'^cm-^)^ which 
may make it more difficult to detect than DCO'*". The 
ground-state line of ortho-H2D"'" at 372 GHz was searched 
toward three sources (LkCa 15, TW Hya and MWC 480) 
in a setting together with N2H+, but neither was detected. 
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Table 5. Calibration sources used for the JCMT observations 



Line 


LkCal5 


TWHya 


HD 163296 


MWC 480 


'"CO J = 3^2 


N2071 IR 


N2071 IR 
IRC+10216 


IRC+10216 
IRAS 16293-2422 


W3(0H) 


"CO J = 3-»2 


W3(OH) 


IRC+10216 


NGC6334I 
G34.3 

IRAS16293-2422 


W3(0H) 


Ci**0 J = 2-.l 


CRL618 








Ci**0 J = 3^2 


W3(OH) 






N1333IR4 


HCO+ J = 4^3 


W3(OH) 


N2071 IR 


G34.3 


W3(0H) 


U^^CO+ J = 4^3 


IRC+10216 


IRC+10216 
N2071 IR 






DCO+ J = 5-^4 


N1333IR4 

IRC+10216 


IRAS16293-2422 






CN J = 3|^2f 


CRL618 


CRL618 


CRL618 

IRAS 16293-2422 


CRL618 
G34.3 


HCN J = 4^3 


CRL618 


IRC+10216 


IRAS16293-2422 


W3(0H) 
CRL618 
OMCl 


HNC J = 4^3 




N2071 IR 






DCN J = 5^4 




N2071 IR 
IRC+10216 






H^^CN J = 4^3 


not available 


IRAS16293-2422 






CS J = 7^6 


not observed 






not observed 


SO J = 88->77 










H2CO J = 3i2->2ii 




N2071 IR 






H2CO J = 5l5->4i4 


CRL618 


N2071 IR 
IRC+10216 


IRAS 16293-2422 


CRL618 


CH3OH J = 7a'^6k 




not observed 






N2H+ J = 4^3 


not available 


not available 




not available 


H2D+ J^lio^lll 


not available 


not available 




not avaialble 



Note.- When a observation is carried out over more than one run, different standard sources are used depending on their 
availability. The entry ... indicates that the source was not observed in that particular line (cf. Table |S|, so no standard was 
needed. 



Because of the poor atmosphere and higher receiver noise 
at this frequency, the limits for both H2D+ and N2H+ 
are not very deep, except toward MWC 480. Recently, 
Ceccarelli et al. (2004) have published the detection of 
the H2D+ 372 GHz hne from the DM Tau disk using the 
Caltech Submillimeter Observatory, together with a tenta- 
tive feature from the TW Hya disk. Their integrated line 
intensity toward TW Hya is / TuBdV = 0.39 ±0.12 K km 
s~^, compared with our 2a limit of 0.20 K km s~^. Taking 
into account the smaller beam dilution in the JCMT beam 
and the measurement uncertainties, the difference between 
these two data sets is about a factor of two. 



4.2. Disks properties and molecular abundances 

4.2.1. Disk mean density 

The mean density can be constrained from line ratios 
of molecules with high dipole moments such as HCO+, 
CN, HCN or H2CO. A simple excitation analysis was 
performed using an escape probability code described in 



Jansen et al. 1)1994111????^ . The code computes the statis- 
tical equilibrium population of the rotational levels given 
the kinetic temperature, volumn density and column den- 
sity. Integrated temperatures of low- J transitions from Qi 
II2001|I were used to complement our high-J data. Both 
sets of data were corrected for beam dilution by multi- 
plying the observed velocity integrated main-beam tem- 
peratures by the size of the beam as listed in Table El 
The analysis for LkCa 15 and TW Hya has been per- 
formed previously by van Zadelhoff et al. (2D(nj using 
HCO"'' and HCN, and takes both the radial and vertical 
density structure of the disk into account. Here H2CO is 
also used as a diagnostic for LkCa 15 adopting the same 
method. Consistent with their results, we find that the 
densities in the regions probed by our observations range 
from 10^ to 10^ cm^"^. This density refers to the re- 
gion where the molecular lines are emitted. The 
fractions of mass in a given density interval for various 
disk models are shown in Fig. 3 of van Zadelhoff et al. 
IPHITII . In all models (Chiang & Goldreich lTMTl D'Alessio 
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Fig. 1. Line profiles observed toward LkCalS. The dashed line indicates the velocity of the source. Note the different 
antenna temperature scales for the different features. 



et al. 119991 Bell et al. I1997|l . most of the gas is located in 10^ cm ■^. Those densities are sufficient for most transi- 
the region of the disk where the density is greater than tions studied here to be thermalized. We refer to the paper 
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Fig. 2. Line profiles observed toward TW Hya. Note the different antenna temperature scales for the different features. 



of van Zadclhoff et al. H2()()l|l for a detailed discussion on 
the disk models, the densities derived from line ratios and 
the disk location where the lines are expected to become 
optically thick. 

4.2.2. Disk mean temperature 

The mean kinetic temperatures are less well constrained; 
the ratio "CO J = 3 ^ 2 / ^^CO J = 1 ^ of 1.35 ± 
0.4 suggests that Tkin -- 20-40 K for LkCa 15 in the re- 
gion where the ^'^CO emission originates, assuming that 



both lines are optically thin (van Zadclhoff et al. I2()()l|l . 
The bulk material where CO emits is therefore on av- 
erage moderately warm and the density is high enough 
that the level populations can be assumed to be thermal- 
ized for most cases. The ratios of 2.4 ± 0.7 for MWC 480 
and 1.7 ± 0.5 for HD 163296 indicate similar tempera- 
ture ranges. From the CO J = 6^5 / J = 3^2 ratios 
presented in Thi et al. H2001|l . it is found that the upper 
layers of disks have temperatures in the range 25-60 K. 
The gas temperatures derived from the H2 data in Thi et 
al. H2001|l are slightly higher for disks around Herbig Ae 
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Fig. 3. Line profiles observed toward MWC 480. Note the different antenna temperature scales for the different features. 



stars than around T Tauri stars, as expected if the disks 
are heated by the radiation from the central star. The 
H2CO J = 3o3 ^ ^02/ J = 322 -^ 221 ratio is potentially a 
good temperature indicator, but the J = 803 -^ 2o2 hne 
has only been detected toward LkCa 15. The limit on the 
line ratio constrains the temperature to be below 200 K. 
A mean temperature of 25 K is adopted in the remaining 
parts of the paper. 



4.2.3. Disk size 



ometer maps, an attempt has been made to infer them 
directly from our model profiles. Two methods have been 
employed. First, since the ^^CO 3-2 emission line is op- 
tically thick, it probes the surface temperature profile of 
the disk (van Zadelhoff et al. l2()()l|l . Using the method de- 
scribed by Dutrey et al. H1997|l an estimate of the disk size 
can be made from the ^^CO J = 3^2 lines: 



Tinb dv = Tox (pSv) 



Rl 



R„ 



D2 



■ cosz 



^l 



(5) 



The disk sizes are important ingredients for comparing where -Rjn and -Rout are the inner and outer radii, 5v is 
the observed column densities with models. Since sizes the local turbulent velocity (between 0.1 and 0.2 km ~^) 
are notoriusly difficult to derive from low S/N interfer- and p a geometrical factor of the order of 1.5. We adopt 
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Fig. 4. Line profiles observed toward HD 163296. Note the different antenna temperature scales for the different 
features. 



here pSv = 0.3 km s'^, i?in=0 AU, Tex=25, 30, and 50 K 
as the mean disk excitation temperature. The values for 
the inclinations i, distances D (in AU) and beam sizes f^a 
are provided in Tabled and [3 The derived disk sizes are 
given in Table [7| Our estimates for Tex=30 K are simi- 
lar to published values except for HD 163296, which we 
find to have twice the size found by Mannings & Sargent 
H1997|l . who measured it directly from their ^^CO 1-0 map. 
Spectra of ^^CO 3-2 emission line were also generated 
using a standard parametric disk model as described by, 
e.g., Beckwith & Sargent (,1993,1 . The code uses a ray- 
tracing method and assumes that the population of the 
rotational levels is in Local Thermodynamic Equilibrium. 
All disks have a power-law density profile of the form 
n{r) — no{r/l AU)^^'^. The exact value of uq cannot be 
constrained by fitting optically thick lines and we assume a 
typical value of 5 10^'^ cm""^ at 1 AU in the mid-plane. The 
disk is in hydrostatic equilibrium in the vertical direction. 
The inner radius is set to 0.01 AU and is not a significant 
parameter. The free parameters of this model are the ex- 
citation temperature Tin at 1 AU, the inclination i and 
the outer radius -Rout ■ For simplicity, the gas temperature 
in the disk is assumed to have a radial profile power in- 
dex of 0.5 and isothermal in the vertical direction. The 
best fits are found using a downhill simplex method (e.g.. 
Press et al. I1997p . Figure |H1 shows the observed spectra 
and their best fits obtained with the parameters reported 
in Table[71 The outer radii found by this ray-tracing model 
are smaller than those from the optical depth model with 
?Gx = 30 K, which can be ascribed to additional contri- 



butions from warmer gas at large radii not taken account 
in the isothermal disk model. Note that Tin and i?out are 
probably degenerate: Table Ogives two sets of values for 
LkCa 15 that can both fit the spectra. Only high signal-to- 
noise spatially resolved interferometer images can lift this 
degeneracy. The larger outer radius (and smaller inner ra- 
dius temperature) is adopted, which is closer to that found 
by direct fitting of interferometric maps (Qi et al. I2003|l . 
The inclinations are consistent with published values (see 
Table El). 

4.3. Column densities 

Given the sizes, molecular column densities can be derived 
from the observed line strengths. Two additional assump- 
tions need to be made: the excitation temperature and 
the line optical depth. The line ratio analysis shows that 
the lines arise from sufficiently high density regions (10^- 
10^ cm~^) that they can be assumed to be thermalized to 
first order, although small deviations are expected in the 
surface layers (see below). Therefore, a single excitation 
temperature of 25 K is adopted to allow easy comparison 
between the disks. 

The optical depth can be estimated from the ratio 
of lines from two isotopologues, assuming that the two 
species have the same excitation temperature. Such data 
are available for a few species and lines, and the results 
arc summarized in Table |H1 It is seen that both the "'^^CO 
and H^^CO+ lines are very optically thick. An alterna- 
tive method is to compare the size of the optically thick 
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Fig. 5. i^CO J = 3^2 maps toward TW Hya, LkCa 15, HD 163296 and MWC 480. 



blackbody which accounts for the line flux to the actual 
disk radius derived from the optically thick ^^CO line. 
We adopt again the approach of Dutrey et al. (1997) and 
rewrite Eq.© as follows: 



i?iinc(AU) = 106.4 (^ 




TTCOSZ 

(6) 
where 9a is the main-beam diameter at half power in arc- 
sec. For a Gaussian shape, the solid angle is given by 
fla = 1.133612. 

Assuming Tox=25 K, all derived radii are significantly 
smaller than the CO disk sizes, except for CN J = 3|— ^2|. 
This would suggest that the lines from species other 
than CO are less thick in the outer parts of the disks 
(R >300 AU) to which our data are most sensitive, but 
we consider the direct determination through the isotopo- 
logue ratios more reliable. 

When the medium is slightly optically thick (r <3) and 
Tgx >> JcMB with TcMB = 2.73 K, the column density 
of the upper level N^ is given by: 



iVu 



Snki^^w fna + n 



(7) 



where v is the frequency of the transition, r is the mean 
optical depth, W — J Tmb dv is the integrated line in- 
tensity expressed in K km s~^, and fla and fig are the 
telescope main-beam and the source solid angles, respec- 
tively. The ratio {ila + fls)/^s is the beam dilution factor. 
Since the inferred disk sizes are much smaller than the 
beam sizes, {fla + fls)/^s — ^a/^s- The Einstein A^i co- 
efficient of the transition in units of s~^ coefficient is given 
by: 



Aui 



(Mtt^v 



4,, 3 



5u 



(8) 



V 3/ic3 

where /i is the dipole moment of the molecule in Debye, 
S is the line strength, and g^ = gjgxgi is the statistical 
weight of the upper level. Finally, the factor 



r^ = 



1 



(9) 



is the escape probability in the so-called Sobolev or Large 
Velocity Gradient approximation. 

The column density in level u is related to the total 
column density N by: 



N„ 



N 



^/rot [^ ex 7 



gj9Kgie 



-E^/T^, 



(10) 
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Table 6. Integrated line intensities 



Line 



/ T„b dv (K km s"^) 



(K) 



(cm"^) (GHz) 



Telescope Beam 

(") 



Cal.° 



LkCal5 TW Hya HD 163296 MWC 480 



i2cO J = 2 — 1 


16.6 


2.7(3 


230.538 


IRAM30m 


10.7 




1.82 








i^CO J = 3-*2 


33.2 


8.4(3 


345.796 


JCMT 


13.7 


yes 


1.17 


1.98 


3.78 


2.88 


i^co J = 3^2 


31.7 


8.4(3 


330.587 


JCMT 


14.3 


yes 


0.39 


0.24 


0.94 


0.57 


C^'^O J = 2^1 


15.8 


2.7(3 


219.560 


JCMT 


21.5 


yes 


<0.20 








C^'^O J = 3-^2 


31.6 


8.4(3 


329.330 


JCMT 


14.3 


yes 


<0.14 








HCO+ J = 4^3 


42.8 


1.8(6 


356.734 


JCMT 


13.2 


yes 


0.26 


1.26 


1.10 


0.35 


h1^CO+ ,7 = 4-^3 


41.6 


1.8(6 


346.998 


JCMT 


13.6 


yes 


<0.13 


0.07 






DCO+ J = 5-^4 


51.8 


3.0(6 


360.169 


JCMT 


13.1 


yes 


<0.10 


0.11 






CN J = 3| — 2| 


32.7 


6.0(6 


340.248 


JCMT 


13.9 


no 


0.67 


1.14 


0.95 


0.29 


HON J = 4^3 


42.5 


8.5(6 


354.506 


JCMT 


13.3 


yes 


0.25 


0.49 


<0.20 


<0.07 


H^^CN ,7 = 4^3 


41.4 


8.5(6 


345.339 


JCMT 


13.6 


no 




<0.04 






HNC 7 = 4^3 


43.5 


8.5(6 


362.630 


JCMT 


13.0 


no 




<0.05 






DCN 7 = 5-^4 


52.1 


4.8(7 


*> 362.046 


JCMT 


13.0 


yes 




<0.03 






CS 7 = 7-»6 


65.8 


2.9(6 


342.883 


JCMT 


14.0 


no 


<0.08 






<0.08 


SO 7 = 8s^77 


87.7 


1.8(6 


344.310 


JCMT 


13.7 


no 




<0.10 






H2CO 7 = 2i2^1ii 


21.9 


1.0(5 


140.839 


IRAM30m 


17.5 




0.17 




<0.10 


<0.40 


H2CO 7 = 3o3^2o2 


21.0 


4.7(5 


218.222 


IRAM30m 


11.3 




0.14 




<0.30 




H2CO 7 = 322^221 


68.1 


2.3(5 


218.475 


IRAM30m 


11.3 




<0.10 






<0.06 


H2CO 7 = 3i2^2ii 


33.5 


4.5(5 


225.697 


IRAM30m 


10.9 




0.10 




<0.30 




H2CO 7 = 3i2^2ii 


33.5 


4.5(5 


225.697 


JCMT 


22.2 


no 




<0.05 






H2CO 7 = 5i5^4i4 


62.5 


1.7(6 


351.768 


JCMT 


13.4 


yes 


0.29 


<0.04 


<0.20 


<0.09 


CH3OH J = 2k->1k 


6.9 


2.6(3 


96.741 


IRAM-30m 


25.4 




<0.05 




<0.03 




CH3OH 7 = 42 ^3i E+ 


45.4 


3.7(4 


218.440 


IRAM30m 


11.3 




<0.10 






<0.20 


CH3OH J = 5k^4:k 


34.8 


4.5(4 


241.791 


IRAM30m 


10.2 




<0.10 




<0.10 




CH3OH J = 7k-*6k 


65.0 


1.3(5 


338.409 


JCMT 


13.9 


yes 




<0.02 






N2H+ 7 = 4^3 


44.7 


4.4(6 


372.672 


JCMT 


12.7 


no 


<0.10 


<0.30 




<0.05 


H2D+ ,7 = lio^lii 


104.3 


1.2(6 


372.421 


JCMT 


12.7 


no 


<0.10 


<0.20 




<0.05 



Note.- The dots indicate not observed. When a hue is not detected, a 2a upper limit on T„ib in a 0.3 km s^^ bin is 
computed and the same profile as the ^'^CO J = 3^2 line is assumed. The beam size (HPBW) is computed for the 
IRAM-30m using the fitting formula flPi3W(")=2460/frcquency(GHz). " Unless specified, the critical densities are 
taken from Jansen 0.995) for Tkin^lOO K assuming optically thin lines. a(b) means a x lO''. ^ Computed using the 
collisional rate coefficients for HCN. '^ Derived assuming the collisional rate coefficients of Peng & Whiteoak ijl993|l . 
Observation of calibration sources before and/or after the object. More details are given in Table O 

Table 7. Disk sizes from the ^^CO J = 3^2 integrated intensities 



Radius (AU) 
Literature" rcx=25 K rex=30 K Tcx=50 K 



Disk model 
i^g.t T{1 AU)" 



LkCa 15 425 

LkCa 15'* 425 

TW Hya 200 

HD 163296 310 

MWC 480 695 



620 


559 


422 


57 


290 


170 








57 


450 


100 


238 


215 


162 


3.5 


165 


140 


778 


701 


530 


65 


680 


180 


722 


650 


492 


30 


400 


170 



"See Table 2 for references 

" Gas temperature at 1 AU. The radius of disk is taken to be 0.01 AU. 

'^ Values adopted in subsequent analysis 

■^ Alternative fit to LkCal5 data 



where <5rot(Tex) is the rotational partition function, gj is 
the rotational statistical weight factor equal to 2 J + 1 for 
diatomic or linear molecules, gx is the _ftr-level degener- 
acy, gi is the reduced nuclear spin degeneracy, and E^ is 
the energy of the upper level expressed in K. For linear 
molecules, gK=gi=^ for all levels. Formaldehyde, H2CO, 
is an asymmetric top molecule with 7=1/2. Combining 



the above equations, the total column density can be ex- 
pressed as: 

1.67 X IQi-^ 



N 



v^i^S 



~ Wrot v-' ex j^ 



Eu/T„, 



1-e- 






■7^ I / Tmb dv 



(11) 
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Fig. 6. Observed (full lines) and simulated (dahes lines) ^^CO J = 3^2 spectra. The observations are normalized to 
the peak values. The simulations are for a single temperature profile disk model (see Table |3 for parameters). 

Table 8. Optical depths estimates from line ratios between isotopologues 





LkCal5 


TWHya 


HD 163296 


MWC 480 


'^CO J = 3^2 


26.8 




8.3 


18.9 


14.9 


"CO J = 3-^2 


0.44 




0.14 


0.31 


0.25 


HCO+ J = 4^3 


19.4 




3.8 






H"CO+ J = 4^3 


< 0.32 




0.06 






HON J = 4^3 


0.75 


< 


5.1 






H^CN J = 4-+3 


0.012 


< 


0.085 







Note.- We assume that the isotopologues have the same excitation temperature and that [^ C]/[^ C]=60. 



This formula is similar to that of Blake et al. (j 19 8 71) and 
Turner l|1991|l but with the introduction of the escape 
probability and beam dilution factor. 

For linear molecules the line strength is equal to the ro- 
tational quantum number J. The rotational energy level 
structure of the two linear molecules CN and HCN are 
more complicated than those for CO. The spin of the un- 
paired electron for CN (S=l/2) and the nuclear spin of 
^''N (1=1) lead to fine- and hyperfine splitting of the ro- 



tational levels. The observed CN 340 GHz line is a blend 
of three lines that account for ~ 55% of all the hyperfine 
lines arising from the level J=3 fSimon ri997|l . The line 
strength is taken from Avery et al. (199^ . An advantage 
of the hyperfine splitting is that it decreases the optical 
depth in each individual component. The HCN J = 4— >3 
line is also a blend of hyperfine lines but we assume that 
all the flux is included in the observed line. Other con- 
stants used to derive the column densities are taken from 
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existing catalogs (Pickett et al. ll998|l and are summarized 
in Thi H2(JU2afl . The rotational partition functions were 
calculated using the formulae for each molecule in Gordy 
& Cook ifroSHl . 

The radical CN and the molecule HCN have different 
critical densities and the HCN J = 4 ^ 3 line may be sub- 
thermally excited in the upper layer, so that the inferred 
A^(CN)/iV(HCN) ratio varies strongly with density. This 
effect, which can be up to a factor of 2 in the CN/HCN 
abundance ratio, has been corrected using the statistical 
equilibrium calculations described above for the inferred 
range of temperatures and densities. It should be noted 
that this correction assumes that the CN and HCN lines 
come from the same location inside the disks, which is 
probably not the case. In disk models, CN peaks more to- 
ward the lower density surface layers than HCN because 
CN is mostly formed by radical reactions and photodis- 
sociation of HCN (Aikawa et al. I2UU2|I . This effect would 
lead to higher CN/HCN abundance ratios than presented 
here. 

Table 1^1 summarizes the beam-averaged column densi- 
ties and upper-limits for the observed molecules, adopting 
the disk sizes derived from the fits to the ^^CO J = 3^2 
spectra using the isothermal disk model (see parameters 
in rightmost columns of Table 01 . A single excitation tem- 
perature Tox = Tkin = 25 K and an optical depth of r = 1 
are assumed for all lines. For optically thick lines with 
T ^ 1, such as those of HCO+ and perhaps HCN, the de- 
rived column densities by these formulae are clearly lower 
limits. Wherever available, column densities derived from 
isotopic data have been used in these cases. 

4.4. Molecular abundances 

It is well known that disk masses derived from CO mea- 
surements are much lower than those obtained from mil- 
limeter continuum observations assuming a gas/dust ratio 
of 100. Because CO is subject to photodissociation and 
freeze-out, one cannot adopt the canonical CO abundance 
of CO/H2=10~'' found for molecular clouds; instead, the 
disk masses Mdisk are assumed to be given by the millime- 
ter continuum observations (see Thi et al. 2001, Table [21). 

It should noted that, in the optically thin limit, the 
abundances are independent of the disk size. The derived 
abundances are summarized in Table [Till As noted above, 
the abundances derived from the highly optically thick 
HCO+ and HCN lines are likely to be underestimated by 
up to an order of magnitude. For molecules that are de- 
tected in all four disks (CN, HCO+), the abundances vary 
significantly from object to object. The non-detection of 
HCN toward the Herbig Ae stars confirms the low abun- 
dances in these cases, although the high critical density of 
the HCN J = 4 ^ 3 line may also play a role. The upper 
limits are much lower in the case of TW Hya owing to the 
small distance of this object and its narrow lines. 

The derived abundances and depletion factors are 
roughly consistent with the large range of values given 



by van Zadelhoff et al. (2001), especially for their colder 
models. A full comparison between the two studies is dif- 
ficult since van Zadelhoff et al. performed more detailed 
radiative transfer modeling with a varying temperature in 
the vertical direction. Also, they adopted a smaller disk 
radius of 200 AU compared with our value of 450 AU for 
LkCa 15. This leads to higher abundances to reproduce 
the same line flux, at least for optically thick lines. 

For comparison, we have also re-derived the abun- 
dances in the disk around DM Tau, where many of the 
same species have been detected by Dutrey et al. H1997|l . 
Their tabulated velocity integrated flux densities have 
been converted to velocity integrated main beam temper- 
atures via the relation: 



J T^b{K) dv = W-^^ X J F{3y) 



dv 



A^(cm) 
2fc 



n- 



(12) 



where F is the flux density in Jansky, k is the Boltzmann 
constant in erg K~^, A the wavelength in cm, Via is the 
main-beam solid angle and dv is in km s^^. A total 
disk mass of 0.018 Mq has been used, computed using 
Equation (6) of Thi et al. H2001|l and a continuum flux 
of 110 mJy at 1.3 mm (Guilloteau & Dutrev lTUMjI . The 
abundances are reported in Table lTUI Our new abundances 
estimates are within a factor of 4 of those deduced by 
Dutrey et al. (|1997fl who used a different method to de- 
rive their abundances. 

The last column of Table ^1 contains the abundances 
found in the cold outer region of the protostellar envelope 
of IRAS 16293-2422. The latter abundances seem to be 
higher than those in disks by at least an order of magni- 
tude, even taking into account the fact that the disk abun- 
dances may be underestimated because of optical depth 
effects. A noticeable exception is CN, which has a higher 
abundance in all four disks and in DM Tau than in IRAS 
16293-2422. 

5. Discussion 

5.1. Molecular depletion 

Two processes have been put forward to explain the low 
molecular abundances in disks. First, in the disk midplane, 
the dust temperature is so low (< 20 K) and the density 
so high (nn > 10^ cm""^) that most molecules including 
CO are frozen onto the grain surfaces. This possibility is 
supported by the detection of large amounts of solid CO at 
infrared wavelengths in the disk around the younger class I 
object CRBR 2422.8-3423 (Thi et al. 2002b). In this en- 
vironment, surface chemistry can occur but the newly- 
formed species stay in the solid phase and thus remain 
unobservable at millimeter wavelengths, except for a small 
fraction which may be removed back in the gas phase by 
non-thermal desorption processes such as cosmic-ray spot 
heating. 

Second, the photodissociation of molecules in the up- 
per layers of protoplanetary disks by the ultraviolet radia- 
tion from the central star and from the ambient interstellar 
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Table 9. Beam averaged column density of various molecules using the disk sizes listed in Tabled 









N{cm-''') 






LkCal5 


TWHya 


HD 163296 


MWC 480 


Species 


Transition 


i?disk=450 AU 


i?disk=165 AU 


i?disk=680 AU 


i?diBk=400 AU 


i^^Qa 


"CO J = 3^2 


1.9(16) 


3.2(16) 


3.5(16) 


6.9(16) 


"CO 


J = 3^2 


3.6(14) 


5.5(14) 


5.9(14) 


1.2(15) 


HCO+ 


J = 4-^3 


3.3(11) 


4.4(12) 


9.4(11) 


1.0(12) 


HCO+ 


H"CO+ J = 4^3 




1.2(13) 






H"CO+ 


J = 4^3 


<1.5(11) 


2.0(11) 






DCO+ 


J = 5^4 


<2.9(11) 


4.4(11) 






CN 


J = 3|-2| 


1.5(13) 


6.6(13) 


1.5(13) 


1.5(13) 


HCN 


J = 4^3 


1.8(12) 


9.2(12) 


<1.0(12) 


<1.2(12) 


H"CN 


J = 4^3 




<4.8(11) 






HNC 


J = 4^3 




<1.4(12) 






DCN 


J = 5^4 




<4.0(10) 






cs 


J = 7->6 


5.1(12) 








Ha CO 


J = 2i2^1ll 


5.1(12) 




<2.1(12) 


<2.7(13) 


H2CO 


J = 3o3 — >2o2 


1.7(12) 




<2.6(12) 




H2CO 


J = 322 ^221 


<1.4(13) 




<9.4(12) 




H2CO 


J = 3i2 -^2ii 


7.1(11) 




<1.5(12) 


<9.4(11) 


H2CO 


J=5l5^4i4 


2.4(12) 


<8.0(11) 


<1.1(12) 


<1.5(12) 


CH3OH 


J^2k^1k 


<7.1(13) 




<2.8(13) 




CH3OH 


J = 42^3i E+ 


<4.3(14) 






<6.9(13) 


CH3OH 


J = 5k^4k 


<2.4(13) 








CH3OH 


J = 7k^6k 




<1.1(13) 






N2H+ 


J = 4^3 


<1.4(12) 


<1.0(13) 




<1.5(12) 


H2D+ 


J=ll0^1ll 


<8.7(11) 


<4.4(12) 




<1.0(12) 


SO 


J = 88^77 




<4.3(12) 







Note.- The excitation temperature is assumed to be 25 K for all lines; uncertainties in the column densities are of the order of 

30%, not including uncertainties in the disk size. 

a(b) means ox 10*". 

" ^^CO column density derived from "CO intensity assuming ["C]/["C]=60. 



medium can limit the lifetime of molecules. The ultravi- 
olet flux from the central star can reach 10^ times the 
interstellar flux (Glassgold, Feigelson & Montmerle I2UUU|I . 
Aikawa et al. H2Ul)2|l and van Zadelhoff et al. H2Ul)3|l have 
modeled the chemistry in disks, taking these mechanisms 
into account. Their models show that molecules are abun- 
dant in the intermediate height regions of disks, consis- 
tent with the derived temperature range (20-40 K) for 
the emitting gas. According to the flaring disk model, this 
intermediate region is located just below the warm upper 
layer (T ~ 100 K). 

The molecular abundance distributions predicted by 
the above chemical models including photodissociation 
and freeze-out have been put into a 2-D radiative trans- 
fer code to compute the level populations using statisti- 
cal equilibrium rather than LTE and to take the optical 
depth effects properly into account. The resulting inte- 
grated fluxes can be compared directly with observations. 
As shown by Aikawa et al. H2002|l they differ by factors of 
a few up to an order of magnitude, which indicates that 
such models are to first order consistent with the data. 



5.2. CN/HCN abundance ratio 

Table ^2 includes the CN/HCN abundance ratios de- 
rived for the disks. Compared with IRAS 16293-2422, the 
CN/HCN ratio is more than two orders of magnitude 
larger, and even compared with galactic PDRs, all disk 
ratios are higher. The disk ratios may be overestimated 
due to underestimate of HCN optical depth effects, but 
the high values are a strong indication that photodisso- 
ciation processes play a role in the upper layers of the 
disks. 

CN is particularly enhanced by photochemistry since 
it is produced by radical reactions involving atomic C 
and N in the upper layers as well as photodissociation 
of HCN. Moreover, CN cannot be easily photodissociated 
itself since very high energy photons (<1000 A, > 12.4 eV) 
are required to destroy the radical (van Dishoeck 11987(1 . 
The CN/HCN ratio appears to be higher in disks around 
Herbig Ae stars than around T Tauri stars, although our 
high ratio of > 11 for MWC 480 disagrees with the ratio 
of -4 by Qi PUUTJI with OVRO. The disagreement be- 
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X = iV/iV(H2) 






Species 


LkCal5 


TWHya 


HD 163296 


MWC 480 


DM Tau 


IRAS 16293-2422'' 




This work" 


Dutrey et al. 




CO 


3.4(-07) 


5.7(-08) 


3.1(-07) 


6.9(-07) 


9.6(-06) 


1.5(-05) 


4.0(-05) 


HCO+ 


5.6(-12) 


2.2(-ll)'= 


7.8(-12) 


I.O(-IO) 


7.4(-10) 


7.4(-10) 


1.4(-09) 


H"CO+ 


<2.6(-12) 


3.6(-13) 






<3.6(-ll) 


<3.6(-ll) 


2.4(-ll) 


DCO+ 


<2.31(-12) 


7.8(-13) 










1.3(-11) 


CN 


2.4(-10) 


1.2(-10) 


1.3(-10) 


1.4(-10) 


9.0(-09) 


3.2(-09) 


8.0(-ll) 


HCN 


3.1(-11) 


1.6(-11) 


<9.1(-12) 


<1.1(-11) 


4.9(-10) 


5.5(-10) 


l.l(-09) 


H"CN 




<8.4(-13) 










1.8(-11) 


HNC 




<2.6(-12) 






1.5(-10) 


2.4(-10) 


6.9(-ll) 


DCN 




<7.1(-14) 










1.3(-11) 


CS 


<8.5(-ll) 








2.4(-10) 


3.3(-10) 


3.0(-09) 


H2CO 


4.1(-11) 


<1.4(-12) 


<1.0(-11) 


<1.4(-11) 


2.4(-10) 


5.0(-10) 


7.0(-10) 


CH3OH 


<3.7(-10) 


<1.9(-11) 


<1.5(-10) 


<2.0(-09) 






3.5(-10) 


N2H+ 


<2.3(-ll) 


<1.8(-11) 




<1.5(-11) 


<5.0(-09) 


<2.0(-10) 




H2D+ 


<1.5(-11) 


<7.8(-12) 




<1.0(-11) 








SO 




<4.1(-11) 










4.4(-09) 



° Re-analysis of data from Dutrey et al. II1997L see text 
^ Outer envelope abundances from Schoier et al. (|2()()2|l 
" Value inferred from 11^^00+ 



Table 11. Relative molecular abundances 



Source 


A/-(HCO+) 


iV(CN) 


iV(DCO+) 




N{CO) 


Af(HCN) 


A/-(HCO+) 


LkCa 15 


1.6(-5) 


7.9 


<0.411 


TWHya 


3.8(-4)i 
1.4(-4)2 


7.1 


0.035 


DM Tan" 


5.3(-5) 


5.8 




DM Tan*" 


7.6(-5) 


18.4 




HD 163296 


2.7(-5) 


>12.4 




MWC 480 


1.5(-5) 


>11.7 




IRAS 16293-2422^= 


3.6(-5) 


0.07 




TMC-T' 


1.0(-4) 


1.5 




Orion Bar'' 


2.0(-5) 


3.8 




IC 63^ 


2.7(-5) 


0.7 





References, (a) Dutrey et al. I|1997|l : (b) from Table 7; (c) Schoier et al. (2002); (d) Ohishi et al. p992|l : (e) Hogerheijde et al. 

lIT^^ : (f) Jansen et al. lIT^^ . 

Notes. 1 Value derived from H"C0+, assuming [^2C]/["C]=60. ^ Value derived from HCO+. 



tween the results of Qi 1)200 l|l and ours can be ascribed to 
the fact that Qi (2001) had H^^cN data available to con- 
strain the optical depth of the HCN line. Also, their 1-0 
lines are less sensitive to the adopted disk density struc- 
ture. In general the fluxes from MWC 480 for transitions 
which have high critical densities are lower than those for 
T Tauri stars, whereas the CO fluxes (with lower critical 
densities) are higher. This may imply that the level pop- 
ulations are subthermal for the disks around Herbig Ae 
stars. 



Van Zadelhoff et al. f2Q{W\ have investigated the effects 
of different UV radiation fields on the disk chemistry, fo- 
cusing on T Tauri stars with and without excess UV emis- 



sion. CN is clearly enhanced in the upper disk layers for 
radiation fields without any excess UV emission owing to 
its reduced photodissociation. When convolved with the 
JCMT beams, however, the differences are difficult to dis- 
cern: the HCN emission is nearly identical for the differ- 
ent radiation fields, whereas the CN emission varies by 
only a factor of a few. Since T Tauri stars like TW Hya 
have been observed to have excess UV emission (Costa 
et al. 12000(1 — probably originating from a hot boundary 
layer between the accretion disk and the star — , the differ- 
ence in the CN/HCN chemistry with the Herbig Ae stars 
may be smaller than thought on the basis of just the stel- 
lar spectra. Bergin et al. H2008|l suggest that strong Lya 
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emission dominates the photodissociation rather than an 
enhanced continuum flirx. Since CN cannot be photodisso- 
ciated by Lya radiation but HCN can (Bergin et al. 2003, 
van Zadelhoff et al. 2003), the CN/HCN ratio is naturally 
enhanced. 

Other chemical factors can also affect the CN/HCN ra- 
tio. Radicals such as CN are mainly destroyed by atomic 
oxygen in the gas-phase and therefore a lower oxygen 
abundance can increase the CN/HCN ratio. Since atomic 
oxygen is a major coolant for the gas, a lower abun- 
dance will also maintain a higher mean kinetic temper- 
ature. Alternatively, the dust temperature could be in 
the regime that HCN is frozen out but CN not because 
the two molecules have very different desorption energies 
(.Edes(CN) = 1510 K and £;dcs(HCN)=4170 K: Aikawa et 
al.[TM3. 

Yet an alternative explanation for high CN abundances 
is production by X-ray photons emitted from the active 
atmosphere of T Tauri stars (e.g., Aikawa & Herbst n"999al 
Lepp & Dalgarno ll996|l . TW Hya is a particularly strong 
X-ray emitter, with a measured X-ray flux 10 times higher 
than the mean X-ray flux observed toward other T Tauri 
stars (Kastner et al.[2002). In addition, H2 w = l->0 S(l), 
another diagnostic line of energetic events, has been ob- 
served toward this object (Weintraub, Kastner, & Bary 
I2000|l . TW Hya may however constitute a special case 
since neither LkCa 15 nor DM Tau seems to show en- 
hanced X-ray emission, yet they have a similar CN/HCN 
ratio. Further observations of molecules in disks around 
strong X-ray emitting pre-main-sequence stars are war- 
ranted to better constrain the contribution of X-rays on 
the chemistry in disks. 

5.3. HCO'^/CO abundance ratio 

Table ^2 compares the HCO+/CO ratios found in the 
disks with those found in a protostellar region (IRAS 
16293-2422), a dark cloud (TMC-1) and two galactic 
photon-dominated regions (PDRs) (Orion Bar and IC 63). 
Within a factor of two, all values are very similar, except 
for the TW Hya disk. It should be noted, however, that 
except for TW Hya, the ratios in disks have been derived 
from the optically thick HCO+ line and may therefore 
be underestimates. Indeed, the ratio obtained using the 
main HCO+ isotope for TW Hya is closer to that of the 
other objects. Observations of H^^CO+ for all disks are 
warranted to make definitive conclusions. 

HCO+ is produced mainly by the gas phase reaction 
Hj^ -|- CO -^ HCO+ -|- H2. Its formation is increased by 
enhanced ionization (e.g., by X-ray ionization to form H3 
in addition to cosmic rays) and by enhanced depletion 
(which also enhances H^ , see e.g., Rawlings et al. |T992). 
The fact that all HCO+ abundances in disks are higher 
than those found in normal clouds (after correction for 
HCO"*" optical depths) suggests that these processes may 
play a large role in the intermediate warm disk layer where 
both molecules are thought to exist. In this context it is 



interesting to note that TW Hya has the largest deple- 
tion of CO and is also the most active X-ray emitter (see 
below). 

The derived HCO+ abundances provide a lower limit 
to the ionization fraction in disks. The typical values of 
10^^^ — 10^^° are high enough for the magnetorotational 
instability to occur and thus provide a source of turbulence 
and mixing in the disk (e.g., Nomura I2002|l . Ceccarelli et 
al. ij2004(l used their H2D+ observations toward DM Tau 
and TW Hya to derive an electron abundance of (4-7). 
10^^", assuming that Hj^ and H2D+ are the most abun- 
dant ions. These values refer to the midplane of those disks 
where the depletion of CO enhances the H2D+/H;^ ratio 
(see § 5.5), whereas our values apply to the intermediate 
layer where most of the HCO+ emission arises. 



5.4. H2CO/CH3OH 

Lines of H2CO have been detected toward only one source, 
LkCa 15. CH3OH is not detected in the single-dish obser- 
vations of any disks, although it is seen in the OVRO 
interferometer data toward LkCa 15 by Qi (2001), who 
derived a CH3OH column density of 7-20 x lO^'' cm^^ 
compared with our upper limit of 9.4 x 10^^ cm^^. Our 
upper limits for H2CO and CH3OH are derived from the 
spectra with the lowest r.m.s. Methanol has been detected 
from the class protostar L1157 by Goldsmith, Langer & 
Velusamy (1999), where its emission has been ascribed to a 
circumstellar disk. However, this object is much younger 
than those studied here and presumably has a different 
physical structure and chemical history. 

The H2CO/CH3OH abundance ratio of > 0.15 for 
LkCa 15 is consistent with values found for embedded 
YSOs (see Table 8 of Schoier et al. IMT for IRAS 
16293-2422 and van der Tak et al. ITOOO for the case 
of massive protostars). Heating of the disk, whether by 
ultraviolet- or X-rays, should lead to strong ice evapora- 
tion and thus to enhanced gas-phase abundances for grain- 
surface products. Both CH3OH and H2CO have been de- 
tected in icy mantles, with the CH3OH abundance vary- 
ing strongly from source to source (Dartois et al. 1999, 
Keane et al. lMITl Pontoppidan et al.|2M|l- For the few 
sources for which both species have been seen, the solid 
H2CO/CH3OH ratio varies from 0.1-1. Thus, the ob- 
served ratio in LkCa 15 could be consistent with grain 
surface formation of both species. Since their absolute 
abundances arc much lower than typical ice mantle abun- 
dances of 10^^ with respect to H2, this indicates that most 
of the CH3OH and H2CO, if present, is frozen onto grains. 
Deeper searches for both species in disks are warranted. 

Protoplanetary disks are places where comets may 
form and their volatile composition may provide con- 
straints on their formation models. The CH3OH abun- 
dance is known to vary significantly between comets. For 
example, comet C/1999 HI (Lee) shows a CO/CH3OH 
ratio around 1 whereas Hale-Bopp and Hyakutake have 
ratios of 10 and 14 respectively (Biver et al. 12000(1 . 
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Comet Lee probably belongs to the so-called "methanol- 
rich comets" group (Bockelee-Morvan, Brooke & Crovisier 
119981 Davies et al. 11993(1 . In addition, the measured CO 
abundance is ~1.8± 0.2% compared to H2O, 5 times less 
than found in Hale-Bopp. Alternatively, Mumma et al. 
H2()()ljl propose that Comet Lee has been heated suffi- 
ciently after its formation for CO to evaporate but not 
CH3OH, so that CH3OH abundance is not enhanced 
but rather CO is depleted. Mumma et al. H^UOlfl notice 
that the CH3OH/H2O and CO/H2O ratios vary strongly 
among comets coming from the giant-planets regions. The 
picture is not complete since CO can be converted to 
CO2, whose abundance is high in interstellar ices (e.g., 
Ehrenfreund & Charnley 20^01) but less well known in 
comets (10% in comet 22P/Kopfr, Crovisier et al.il999j). 

Long period comets were probably formed in the 
Jupiter-Saturn region (around 5-20 AU), whereas our 
data are only sensitive to distances of more than 50 AU. 
It would therefore be more relevant to compare the com- 
position of protoplanetary disks to that of Kuiper Belt 
Objects, which were formed beyond 50 AU in the solar 
nebula. The chemical composition of Kuiper Belt Objects 
is not well known (see Jewitt & Luu l2000|l , although obser- 
vations show that comet nuclei and Kuiper Belt Objects 
have different surface compositions (Luu & Jewitt 2003 
Jewitt l?()()2|l . The nature of Centaur objects is better un- 
derstood. It is believed that Centaur objects were formed 
beyond 50 AU and recently entered the planetary zone 
with orbits crossing those of the outer planets. The best 
studied Centaur object, 5145 Pholus, shows the presence 
of CH3OH although the exact amount is not well con- 
strained (Cruikshank et al. I1998|l . 

5.5. D/H ratio in a circumstellar disk 

The D/H ratio in the TW Hya disk of 0.035 ±0.01 has 
been derived from the H^'^CO"'" and DCO"*" column density 
ratios, assuming an isotopic ratio [-^^C]/[^'^C] of 60 (van 
Dishoeck et al lM)^ . Hints of H"CN and DON features 
are seen in the TW Hya spectra, but neither of them is 
definitely detected. Searches for other deuterated species 
in the LkCa 15 disk, in particular DCN and HDO, are 
reported by Kessler et al. I|2()03|l . 

It should be kept in mind that our observations provide 
only a value of the D/H ratio averaged over the entire disk. 
Models of the DCO"'"/HCO'*' abundance ratio show that 
it remains rather constant throughout the disk down to 
a radius of 100 AU, however (Aikawa et al. ;2002). AU 
values are significantly higher than the elemental [D]/[H] 
abundance of 1.5x10^^ (Pettini & Bowen i2001i Moos et 
aLHHHl. 

Theoretically, the amount of deuterium fractionation 
in molecules depends on the gas kinetic temperature, 
which drives the isotopic exchange reactions, and on the 
cosmic ray ionization rate (Aikawa & Herbst [T999b|l . Also, 
the abundance is enhanced if CO is significantly depleted 
onto grains (Brown & Millar 11989(1 . Thus, the amount 



of deuterium fractionation can serve as a tracer of the 
temperature history of the gas. The deuterium fractiona- 
tion can be further enhanced by grain-surface formation 
(Tielens ,1983j , although not for DCO+/HCO+. Recent 
chemical models succeed in explaining the high fraction- 
ation observed here and in dark cloud cores (Rodgers & 
Millar 1996: Roberts & Millar IMmH Tine et al. lMin|l . but 
only if significant freeze-out is included (Roberts et al. 
I2002L t2003 1. Our observed values are also close to those 
found in disk models which include a realistic 2D temper- 
ature and density profile with freeze-out (Aikawa et al. 

mm- 

Table 2 in van Dishoeck et al. (2003) compares the 
D/H ratio found in disks to typical values for the D/H 
ratio in different protostellar and cometary environments. 
The value found in disks is somewhat higher than that 
in the low-mass protostellar envelope of IRAS 16293- 
2422, but comparable to that seen in dark cloud cores. 
DCO"'"/HCO"'" has not been observed in comets, but the 
D/H ratios derived from DCN/HCN in pristine mate- 
rial in jets originating from below the comet surface is 
found to be similar to that seen for DCO+/HCO+ in the 
TW Hya disk (Blake et al. 1999). Alternatively, Rodgers 
& Charnley (,2002,1 propose that the DCN and HON seen 
in these cometary jets are the photodestruction products 
of large organic molecules or dust grains. In either case, 
the D/H ratio of pristine icy material in comets is high. 
The similarity suggests that either the gas is kept at low 
temperatures as it is transported from pre-stellar clouds 
to disks and eventually to comet-forming regions, or, al- 
ternatively, that the DCO"''/HCO+ ratio has been reset 
by the chemistry in disks. Comparison with D/H ratios of 
molecules which likely enter the disks as ices are needed 
to distinguish these scenarios. 

6. Conclusions 

We surveyed low- and high- J transitions of simple organic 
molecules in two classical T Tauri and two Herbig Ae stars. 
Analysis of line ratios indicates that the emission comes 
from a dense (nn = 10^-10* cm~'^) and moderately warm 
region (T ~ 20-40 K). Detailed fits to the ^^CO 3-2 emis- 
sion line profiles provide independent estimates of the sizes 
of the disks. 

Emission from the ion HCO^ and the radical CN are 
particularly strong, indicating an active gas-phase chem- 
istry in the surface layers of disks which is affected by UV 
radiation from the central stars. H2CO is detected in one 
source but CH3OH is not observed in any object in our 
sample. In one source (TW Hya) the detection of DCO+ 
allows to constrain the DCO^/HCO"*" ratio to ~0.035, a 
value that is higher than that found in the envelopes of 
low-mass protostars but comparable to that observed in 
cold dark cores, where fractionation due to low tempera- 
ture chemistry and CO freeze-out is important. 

This work demonstrates that organic chemistry in 
disks around low- and intermediate-mass pre-main- 
sequence stars can now be studied observationally. The 
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detection of molecular species in disks is hampered, how- 
ever, by the small sizes of disks compared with the actual 
beams of single-dish telescopes. Moreover, because the to- 
tal amount of material is small (few x 10~^ Mq), the 
observations are limited to the most abundant species. It 
is likely that the chemistry in more tenuous disks, in which 
the ultraviolet radiation can penetrate through the entire 
disk, is different from that for our objects (e.g., Kamp & 
BertoldiEOnOl Kamp et al. PTO^ . Although the outer disks 



can be resolved by current millimeter interferometers, in- 
tegration times are too long to do molecular line surveys 
and the inner tens of AU are still out of reach. The detec- 
tion of more complex and much less abundant molecules 
in protoplanetary disks at different stages of evolution 
awaits the availability of the Atacama Large Millimeter 
Array (ALMA). Complementary infrared observations of 
solid-state species along the line of sight of edge-on pro- 
toplanetary disks will help to constrain quantitatively the 
level of depletion in the mid-plane of disks. 
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